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ABSTRACT 

Using Greenwich and SOON sunspot group data during the period 1874-2005, we 
find that the sums of the areas of the sunspot groups in 0° - 10° latitude-interval of 
the Sun's northern hemisphere and in the time-interval, minus 1.35 year to plus 2.15 
year from the time of the preceding minimum-and in the same latitude interval of 
the southern hemisphere but plus 1.0 year to plus 1.75 year from the time of the 
maximum-of a sunspot cycle are well correlating with the amplitude (maximum of 
the smoothed monthly sunspot number) of its immediate following cycle. Using this 
relationship it is possible to predict the amplitude of a sunspot cycle by about 9-13 
years in advance. We predicted 74 ± 10 for the amplitude of the upcoming cycle 24. 
Variations in solar meridional flows during solar cycles and 9-16 year variations in 
solar equatorial rotation may be responsible for the aforementioned relationship. 
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meridional flows (jZaatri et al . 2006). Therefore, north-south 
asymmetry in solar activity is an important physical solar 
property and it greatl y helps for understanding variations 
in the solar activity ([Sokolofj 1 19941 : Ijavaraiah fc Gokhald 
Il997al : iKnaack et al.ll2005l ). In this letter we have used this 
property of a solar cycle to predict the amplitude of the 
upcoming solar cycle 24. 



2 DATA ANALYSIS AND RESULTS 

We have used the Greenwich sunspot group data during the 
period 1874-1976, and the sunspot group data from the 
Solar Optical Observing Network (SOON) of the US Air 
Force/US National Oceanic and Atmospheric Administra- 
tion during 1977 January 1 - 2005 September 30. We have 
taken recently updated these data from the NASA web- 
site of David Hathaway (http : / / solarscience .msf c .nasa. 
gov/greenwich. shtml). These data include the observation 
time (the date and the fraction of the day), the helio- 
graphic latitude and the longitude, central meridian distance 
(CMD), the corrected whole spot area (in mh), etc. for each 
day of the spot group observation (130 mh w 10 22 Mx). 
In the present analysis we have excluded the data corre- 
sponding to the \CMD\ > 75° in any day of the spot group 
life-time. This precaution considerably reduces the errors 
in the derived results due to the foreshortening effect. In 
case of SOON data, we increase area by a factor of 1.4. 
David Hathaway found this correction is necessary to have 
a combined homogeneous Greenwich and SOON data (see 



1 INTRODUCTION 

The prediction of the level of activity is i mportant because 
solar activity impact us in ma ny ways l|Hathawav et. all 
Il999l : IHathawav fc WilsonI l2004n . For example, solar flare 
activity cause geomagnetic storm that can cripple commu- 
nication and damage power grids. There is also mounting 
evidence that solar activit y has an influence on terrestrial cli- 
mate and space weather llRozelot| [2001; Hir emath fc Mandil 
|2004 iGeorgieva et~aT1 120051 ). Many attempts have been 
made to predict the amplitude of a new sunspot cycle by 
using old cycles data with a bel ief that solar magneti c 
field persists for quite sometime (|Hathawav et. all fl999). 
The existence of a statistically significant difference be- 
tween the levels of solar activity in the northern and the 
southern hemispheres is shown by several st atistical stud- 
ies for most of the so lar activity phenomena l|Garcialll990l ; 
ICarbonel etal]|l993h . The north-south asy mmetry is un - 
usually large during the Maunder minimum l|Sokolofi1ll994l ). 
The existence of a few periodicities in the north-south asym- 
metry of solar activity is a lso shown (Ijavaraiah fc Go khalc 
Il997al : IKnaack et al.ll2005l ). In addition, there are consider- 
able north-south differences in the differ ential rotation rates 
and t he meridional motions of sunspots l|javaraiah fc Ulrichl 
2006). Helioseismology measurements also show the exis- 
tence of north-south differences in the solar rotational and 
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aforementioned web-site of David Hathaway.) We binned the 
daily data into 10° latitude intervals, in both the northern 
and the southern hemispheres, and determined sum of the 
areas (AT) of the spot groups in each 10° latitude inter- 
val, separately for the rising and the declining phases of the 
sunspot cycles 11 - 23. (It should be noted here that cycle 23 
is not yet complete. The data are available for about 9 years 
of this cycle. In case of cycle 11, the data are available for 
the last 4 years.) 

We determined cross-correlations between AT and 
amplitude of cycle (RM). We have taken the values of 
RM (which is the largest smoothed monthly mean sunspot 
number), and the epochs of maxima (TM) and the pre- 
ceding minima (Tm) of cycles 12-23 from the web-site, 
|f tp : lit tp . ngdc . noaa . gov/STP/SOLAR_DATA /SUMSPOT_HUMB 
ERS. Fig. 1 shows the cross-correlation function, 
CCF(RM, AT), in different latitude intervals (a posi- 
tive value of lag indicates that RM leads AT). In this figure 
it can be seen that except for AT during the declining 
phases of the cycles and in 0° - 10° latitudes intervals of 
both the northern and the southern hemispheres, for each 
of the remaining cases, viz. AT in 10° -20° and 20° -30° 
latitude intervals during the declining phases of the cycles 
and in all the latitude intervals during the rising phases of 
the cycles, the corresponding CCF(RM, AT) has a weak 
peak at lag ^ 0. This suggests that in all these latitude 
intervals AT and RM variations are approximately in the 
same phase or RM leads AT. In case of AT during the 
declining phases of the cycles, in 0° - 10° latitude interval 
of the southern hemisphere the CCF(RM, AT) has a well 
defined peak (value 0.76) at lag — —1, suggesting that AT 
leads that of RM by about 5-10 years. In the same latitude 
interval of the northern hemisphere the CCF(RM, AT) is 
found to be having a broad peak with two humps (values 
0.8 and 0.6) at lag = and lag — —2, suggesting that AT 
leads RM by about 5-25 years. These results indicate that 
AT can be used to predict RM. 

There exist a number of short-term periodicities, a few 
days to a few years, in bo t h the solar activity and the solar 
rotation (iBai fc Sturrockl 1 19931 ; IJavaraiah fe Komml 1 19991 ; 
iKnaack et al.l l2005). Amplitude of such a periodicity largely 
varies during a solar cycle. Therefore, there is a possibil- 
ity that AT in 0° - 10° latitude intervals of the northern 
and the southern hemispheres during some short intervals 
having strong correlations with RM. With this hypothesis 
we determined the maximal values of correlations between 
AT of cycle n and RM of cycle n + 1 in the following way, 
where n — 12, ...,22 is the cycle number: First we deter- 
mined the values of AT in the intervals which were chosen 
arbitrarily around the epochs of the maxima and the pre- 
ceding minima of the cycles. The AT determinations are 
repeated by increasing or decreasing the lengths of the in- 
tervals with a step of ^ 0.05 year at a time. We find that 
in 0° - 10° latitude interval of the southern hemisphere, the 
correlation is maximum, coefficient of correlation r = 0.97 
(from eleven data points), in the short (0.75 year) time- 
interval just after 1-year after the time of maximum of each 
of the cycles 12-23, TM* : TM + (1.0 to 1.75) (i.e., close 
to the time of the reversal of polarities of the polar mag- 
netic fields). We also find that in 0° - 10° latitude interval of 
the northern hemisphere r = 0.95 is maximum in the time- 
interval (3.5 year), Tm* : Tm+ (-1.35 to 2.15). Both these 



correlations are statistically high significant with > 99.99 
confidence level (from Student's t-test), i.e., the chance of 
getting these relations from uncorrelated quantities is less 
than 0.01%. Interestingly, the e xistence of 0.75 year period- 
icity is known in solar activity (|Knaack et al.ll2005l ). and it 
may be a subharmonic of the well-known Rieg er periodic- 
ity in solar flare activity (|Bai fc Sturrockl |l9 93). The exis- 
tence of 3.5 year periodicity in solar activity is also known 
and this periodicity seems to be more pronounced in the 
north-s outh asymmetries of solar activity and surface ro- 
tation |javaraiah fc Gokhaldll997al ; lKnaack et alj|2005h . In 
Table 1 we have given the values of AT during Tm* and 
TM* . In the same table we have also given the values of the 
amplitudes and the epochs of maxima and minima of the 
sunspot cycles 12-23. 

We find the following linear regressions fits between AT 
and RM correspond to the correlations above: 

RM n+ i = (1.72 ± 0.19) x AT n (Tm*) + (74.0 ± 7.0), (1) 

RM n+1 = (1.55 ±0.14) x AT n (TM*) + (21.8 ± 9.6), (2) 

where uncertainties in the coefficients are the formal l-o 
(standard deviation) errors from the fit. In equations (1) 
and (2) the slopes are on 9a and llcr levels, respectively. 
That is, they are statistically high significant. Therefore, 
the relationship between AT n and RM n +i is well described 
by these linear equations. It should be noted here that al- 
ways Tm* is associated with 0° - 10° latitude interval in the 
northern hemisphere, whereas TM* is associated with 0° - 
10° latitude interval of the southern hemisphere (for other 
combinations, i.e., TM* with 0°-10° interval of the north- 
ern hemisphere and Tm* with 0° - 10° interval of the south- 
ern hemisphere the values of r found to be mere 0.11 and 
-0.24, respectively). 

Using equations (1) and (2) the amplitudes of the up- 
coming sunspot cycles can be predicted by about 13 years 
and 9 years in advance, respectively. The results of the least- 
square fits are shown in Fig. 2(a). Fig. 2(b) shows the corre- 
lation between the simulated amplitudes (PM) [simulated 
using equations (1) and (2)] and the observed amplitudes 
(RM) of the cycles 13-23. The correlations between PM 
and RM and their levels of significance are the same as those 
of AT n and RM n +i. 

Using equation (1) and (2) we obtained the values 
112 ± 13 and 74 ± 10, respectively, for RM of the upcoming 
cycle 24 (the uncertainty is lo value). The latter is more sta- 
tistically significant than the former. Hence, by using equa- 
tion (2) the amplitude of a cycle can be predicted accu- 
rately by 9 years advance. Therefore, we predict 74 ± 10 
for RM of cycle 2 4 . This is equal to the value predicted 
by ISvalgaard et all l|2005h (see Section 3). The pattern of 
the mean cycle-to-cycle variation of the simulated ampli- 
tudes (PM) obtained using equations (1) and (2) is found 
to be slightly more strikingly resemble with that of RM 
(r = 0.97). From this we get 93±10 for RM of cycle 24. How- 
ever, the difference between the values obtained from equa- 
tions (1) and (2) for cycle 24 is significantly large. The mean 
deviation is at 2<r level. Hence, we do not suggest the mean 
value for RM of cycle 24. Moreover, from equations (1) and 
(2) we can get RM n+1 » 2.1 X AT n (TM* ) -0.6 x AT n (Tm*). 
[This may be a more appropriate representation, because 
this is included both terms, A n (Tm*) and AT n (TM*).\ 
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Figure 1. Plots of the CCF(RM, AT) verses lag during the rising and declining phases of solar cycles 12-13. A positive value of lag 
indicates that RM leads AT. The filled circle-solid curve, triangle- dotted curve, square-dashed curve , and open circle- dash- dotted curve 
represent CCF(RM, AT) in latitude intervals 0° - 10°, 10° - 20°, 20° - 30° and in whole disk, respectively. 



Table 1. The maximum (RM) and the minimum (Rm) amplitudes (the largest and the smallest smoothed monthly mean sunspot 
numbers) of the solar-cycles 12 - 23 and the sum of the areas of spot groups (AT, normalized by 1000) in the intervals Tm* = Tm + 
(—1.35 to 2.15) and TM* = TM + (1.0 to 1.75), where TM and Tm represent the maximum and the preceding minimum epochs of the 
solar cycles, respectively. 



Cycle Minimum Maximum Latitude Int.: 0° - 10° (north) Latitude Int.: 0° — 10° (south) 

n Tm Rm TM RM Tm* AT TM* AT 



12 


1878.9 


2.2 


1883.9 


74.6 


1877.55- 


1881.05 


9.47 


1884.90- 


1885.65 


42.11 


13 


1889.6 


5.0 


1894.1 


87.9 


1888.25- 


1891.75 


3.22 


1895.10- 


1895.85 


32.64 


14 


1901.7 


2.6 


1907.0 


64.2 


1900.35- 


1903.85 


12.98 


1908.00- 


1908.75 


54.64 


15 


1913.6 


1.5 


1917.6 


105.4 


1912.25- 


1915.75 


3.74 


1918.60- 


1919.35 


34.58 


16 


1923.6 


5.6 


1928.4 


78.1 


1922.25- 


1925.75 


33.96 


1929.40- 


1930.15 


75.96 


17 


1933.8 


3.4 


1937.4 


119.2 


1932.45- 


1935.95 


29.96 


1938.40- 


1939.15 


82.01 


18 


1944.2 


7.7 


1947.5 


151.8 


1942.85- 


1946.35 


69.35 


1948.50- 


1949.25 


119.65 


19 


1954.3 


3.4 


1957.9 


201.3 


1952.95- 


1956.45 


15.23 


1958.90- 


1959.65 


53.01 


20 


1964.9 


9.6 


1968.9 


110.6 


1963.55- 


1967.05 


50.31 


1969.90- 


1970.65 


78.28 


21 


1976.5 


12.2 


1979.9 


164.5 


1975.15- 


1978.65 


60.05 


1980.90- 


1981.65 


83.53 


22 


1986.8 


12.3 


1989.6 


158.5 


1985.45- 


1988.95 


29.85 


1990.60- 


1991.35 


67.48 


23 a 


1996.4 


8.0 


2000.3 


120.8 


1995.05- 


1998.55 


21.99 


2001.30- 


2002.05 


33.58 



a indicates the incompleteness of the current cycle 23. 
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Figure 2. Plots of the correlations (a) between the AT (for the values given in Table 1) during the intervals Tm* and TM* correspond 
to cycle n and RM of the cycle n + 1, and (b) between RM and the simulated amplitude PM of cycle n + 1, where n = 12, 22, is the 
cycle number. The straight lines represent the corresponding linear relationships. The values of the correlation coefficient (r) are also 
given. The filled circle and the solid line correspond to the AT during Tm* and the open circle and dotted line correspond to the AT 
during TM* . The cross and triangle represent the values for RM of cycle 24 obtained using AT during Tm* and TM* , respectively, 
and the square represents the corresponding mean value. We predict the value represented by the triangle for RM of cycle 24. 



From this relation we get a much smaller value, 57 ± 13, for 
the amplitude of cycle 24 (r = 0.95). It is somewhat closer to 
the value obtained from equation (2). [The negative sign of 
the coefficient of AT n (Tm*) in the aforementioned relation 
can be attributed to the opposite polarities of the magnetic 
fields at Tm* and TM n +i (in sunspot latitude belt).] 

Each of the above derived values for the amplitude of 
the upcoming cycle 24 is less than the RM of cycle 23. This 
is consistent with the indication that the level of activity 
is n ow at the declining ph ase of the current Gleissberg cy- 
cle (| Javaraiah et al.ll2005l ). From equations (1) and (2) we 
can also get AT n (TM*) w l.UxAT n (Tm*)+33.6. [r = 0.94, 
between the simulated and the observed AT(TM*). Note: 
the residual is quite large in case of cycle 23.] Hence, the 
magnetic field at Tm* may contribute to the field at TM n +\ 
both directly and through influencing the field at TM*. 
There is also a suggestion that when AT n (Tm*) is zero 
the AT n (TM*) is not always zero. This might have hap- 
pened during the late Maunder minimum, when sunspot ac- 
tivity is somewhat more pronounced in the southern hem i- 
sphere than in the northern hemisphere (see ISokolofill 19941 ) . 
[The current cycle 23 will be ending soon. So, using equa- 
tion (1), or using the aforementioned relationship between 
AT n {Tm*) and AT„(TM*) and equation (2), an approxi- 
mate prediction can be made for the amplitude of cycle 25 
in a 3 years time.] 



3 DISCUSSION 

The strength of the preceding minimum is used to predict 
the strength of the maximum of the same cycle. However, 
it seems this methods works better after 1-2 year after the 
start of the cycle, i.e., an accurate prediction is possible 
only by about 3-4 years advance. The same is also true for 
the predict ions based on geomag netic indices as precursor 
indicators (Hat haway et. aUll999T ). 

The magnetic fields at the Sun's pol ar regions are im- 
porta nt ingredient for a dynamo model l|Ulrich fc Bovdenl 
120051 ). The polar field is maximum near sunspot mini- 



mum. Scat ten et alj (1 19781 ) have used, for the first time, the 
strength of the polar fields at the preceding minimum of 
a cycle as a precursor indi cator to the strength of the fol- 
lowing maximum. Recentlv. lSvalgaard et all (|2005l ) analyzed 
the polar fields data during the recent four solar cycles and 
predicted a small amplitude, 75 ± 8, for the upcoming cy- 
cle 24. Obviously from this method the prediction can be 
made only by about 5 years in advance. This method seems 
to be more uncertain and could fail if used too early before 
the start of the cycle l|Svalgaard et alj|2005h . 

iDikpati et ail l|2006t) , by simulating the surface mag- 
netic flux using the guidelines of a dynamo model, predicted 
a large amplitude, 150-180, for cycle 24, i .e., a contradic- 
tion t o the aforementioned prediction by ISvalgaard et al.l 
(2005). This discrepancy implies that the dynamo processes 
are not yet fully unders tood, making prediction more diffi- 
cult (|Tobias et al.ll2006l ). 

Using the well known Gnevyshev-Ohl rule or G-0 
rule l|Gnevvshev fc Ohllll94sl) it is possi ble to predic t only 
the amplitude of an odd numbered cycle (Wilson 1988). This 
is also not always possible because occasionally (for exam- 
ple, recently by the cycles' pair 22,23) the G-0 rule is vio- 
lated. A major advantage of the AT n - RM n +i relationships 
above is that using these the amplitudes of both odd and 
even numbered cycles can be predicted. In addition, this 
new method seems to have a solid physical basis. Interest- 
ingly, the TM* is very close to the e poch when the polar - 
fields polarities reversals take place (|Makarov et alj 120031 ) 
and Tm* is close to the epoch when the magnetic fields 
polarities reversals take place close to the equator, i.e., at 
the beginning of a cycle and conti nuing through the years 
of minimum (iMakarov et al.l|200ll ). This suggest that the 
AT n -RM n +i relationships are related to the 22-year so- 
lar magnetic cycle. It should be noted here that although 
sunspot activity is confined to middle and low latitudes, it 
ma y be caused by the global modes of solar magnetic cy- 
cle (|Gokhale et al.lll992l"|juckettll2003l ). 

Reconnection of the magnetic fields of opposite po- 
larities is believed to be the basic mechanism of flare ac- 
tivity. During Tm* the magnetic field structure seem to 
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be largely quadrupole nature, wh ich is proba bly favor- 
able for X-class flares production (|Garcial ll99Cj ). The so- 
lar meridional flows transport angular momentum and mag- 
netic field from pole to equator and vise-versa, in the con- 
vection zone. The motions of spot groups mimic the mo- 
tions in the convection zone (jjavaraiah fe Gokhald Il997bl : 
iJavaraiah fe Kommlll999h . The mean meridional motion of 
sunspot groups is changing from pole-ward to equator-ward 
rapidly in 0° - 10° latitude interval of the northern hemi- 
sphere and gradually in the same latitude interval of south- 
ern hemi sphere during Tm* and TM* , respectively (see 
Fig. 2 in IJavaraiah fe UlrichlboOfJ ). These results indicate 
a participation of the meridional flows in the magnetic re- 
connection process and the reversals of the polarities of 
magnetic fields during Tm* and TM* . The interceptions of 
the pole-ward and the equator-ward meridional flows may 
be responsible for the quadrupole nature of magnetic fields 
during Tm* . It seems that during rising phases of the cy- 
cles the flare activity is strong in the northern hemisphere 
and weak in the southern hemisphere, and t his is opposit e 
during the declining phases of the cycles (|Garcial j 1990) . 
During the rising phases of the cycles the mean merid- 
ional velocity of spot groups is equator-ward in the north- 
ern hemisphere and pole-ward in the southern hemisphere. 
During the declining phases of the cycles the velocity is 
pole-ward in both hemispheres, but the variation is steep 
in the southern hemisphere, ma inly in 20° - 30° latitude in- 
terval (IJavaraiah fe Ulrichll2006T ). In view of the above infer- 
ences, the north-south asymmetry in solar flare activity may 
be related to the north-south asymmetry in the meridional 
flows. The corresponding losses in the magnetic flux in the 
northern and the southern hemispheres caused by the re- 
connection processes may have a contribution for the north- 
south asymmetries in solar magnetic field and in sunspot- 
activity. 

The lengths of the intervals from the beginnings of 
Tm* and TM* of a preceding cycle to TM of its follow- 
ing cycle vary 14-19 years and 7-11 years, respectively. 
The corresponding mean values are found to be 16 years 
and 9.6 years, respectively. Similar periodicities exist in 
both the equatorial rotation rate and the latitude gradi- 
ent term of t he solar rotation determined from the sunspot 
group data (IJavaraiah &s Gokhald Il997al ; IJavaraiah! [2005; 
iGeorgieva et al.l 120051 ). Therefore, variations in the solar 
meridional flows during solar cycles and 9-16 year varia- 
tions in the solar equatorial rotation may be responsible for 
the AT n - RM n +i relationships above. 



4 CONCLUSIONS 

Using Greenwich and SOON sunspot group data during the 
period 1874-2005 we find that: 

(i) The sum of the areas (AT) of the spot groups in 0° - 
10° latitude interval of the Sun's northern hemisphere dur- 
ing the interval Tm* : Tm+ ( — 1.35 to 2.15) in a cycle is well 
correlated with the amplitude (RM) of its following cycle, 
where Tm is the time (in years) of the preceding minimum 
of the preceding cycle, 

(ii) The AT of the spot groups in 0° - 10° latitude interval 
of the southern hemisphere during the interval TM* : TM + 
(1.0 to 1.75) in a cycle is also well correlated with RM of 



its following cycle, where TM is the time (in years) of the 
maximum of the preceding cycle. 

(iii) Using '(i)' and '(h)' it is possible to predict RM of a 
cycle by about 13 years and 9 years advance, respectively. 

(iv) We predicted 74 ± 10 for RM of cycle 24. 

(v) Variations in solar meridional flows during solar cycles 
and 9-16 year variations in solar equatorial rotation may 
be responsible for the relations '(i)' and '(ii)'. 



ACKNOWLEDGMENTS 

I thank the referee, Dr. Leif Svalgaard, for helpful comments 
and suggestions, and Dr. David H. Hathaway for valuable 
information on the data. I also thank Professor Roger K. 
Ulrich and Dr. Luca Bertello for fruitful discussion, and I 
acknowledge the funding by NSF grant ATM-0236682. 



REFERENCES 

Bai T., Sturrock P. A., 1993, ApJ, 409, 676 
Carbonell M., Oliver R., Ballester J. L., 1993, A&A, 274, 
497 

Dikpati M., de Toma G., Oilman P. A., 2006, Geophys. 

Res. Lett., 33, L05102 
Garcia H. A., 1990, Sol. Phys., 127, 185 
Georgieva K., Kirov B., Javaraiah J., Krasteva R., 2005, 

Planet. Space Sci. 53, 197 
Gnevyshev M. N., Ohl A. I., 1948, Astron. Zh. 25, 18 
Gokhale M. H., Javaraiah J., Kutty K. N., Varghese B. A., 

1992, Sol. Phys., 138, 35 
Hathaway D. H., Wilson R. M., Reichmann E. J., 1999, J. 

Geophys. Res., 104, 22,375 
Hathaway D. H., Wilson R. M., 2004, Sol. Phys., 224, 5 
Hiremath K. M., Mandi P. I., 2004, NewA, 9, 651 
Javaraiah J., 2005, MNRAS, 362, 579 
Javaraiah J., Gokhale M. H., 1997a, Sol. Phys., 170, 389 
Javaraiah J., Gokhale M. H., 1997b, A&A, 327, 795 
Javaraiah J., Komm R. W., 1999, Sol. Phys., 184, 41 
Javaraiah J., Ulrich R. K., 2006, Sol. Phys., 237, 245 
Javaraiah J., Bertello L., Ulrich R., 2005, Sol. Phys., 232, 

25 

Juckett D., 2003, A&A, 399, 731 

Knaack R., Stenflo J. O, Berdyugina S. V., 2005, A&A, 
438, 1067 

Makarov V. I., Tlatov A. G., Sivaraman K. R., 2001, Sol. 
Phys., 202, 11 

Makarov V. I., Tlatov A. G., Sivaraman K. R., 2003, Sol. 
Phys., 214, 41 

Schatten K. H., Sherrer P. H., Svalgaard L., Wilcox J. M., 

1978, Geophys. Res. Lett., 5, 605 
Sokoloff D., Nesme-Ribes E., 1994, A&A, 288, 293 
Svalgaard L., Cliver E. W., Kamide Y., 2001, Geophys. 

Res. Lett., 32, L01104 
Rozelot J. P., 2001, J. Atmo. phys., 63, 375 
Tobias S., Hughes D., Weiss N., 2006, Nature, 442, 26 
Ulrich R. K., Boyden J. E., 2005, ApJ, 620, L123 
Wilson R. M., 1988, Sol. Phys., 117, 269 
Zaatri A., Komm R., Gonzalez H. I., Howe R., Corbard T., 

2006, Sol. Phys., 236, 227 



6 J. Javaraiah 

This paper has been typeset from a Tj^X/ P^TfjX file prepared 
by the author. 



